Transport theory: Spatial coherence of random laser emission 
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Recently random laser reached the stage of technological applicability. They have already been engineered 
as coherent microscope light sources in combination with light transport based disordered lenses. The big issue 
for all kinds of applications is the degree of coherence of the emitted radiation. The lasing spot sizes in different 
regimes may provide different degrees of spatial and temporal coherence and as a consequence they can be 
perfectly tunable light sources for the case that the modal behavior can be controlled easily. In this letter we 
investigate the spatial coherence lengths of different random laser samples theoretically. The samples only vary 
in their filling with spherical ZnO Mie scatterers. Beyond we show, that the scattering mean free paths of random 
lasers are not only a material characteristics and dependent to the filling, instead the mean free paths change in 
depth of the sample and therefor depend on the nonlinear self-consistent gain of the random lasing principle. 

PACS numbers: 



INTRODUCTION 



Random lasers have been investigated now for almost 15 
years, and one would think that they are more or less 'over' 
now. But everytime they are said to be dead a whole bunch of 
new ideas and experiments emerge just like the lasing modes 
emerge on another position in the random laser sample. Both 
aspects for sure guarantee the fascination of random lasers as 
playground for theorists as well as experimentalists, however 
beyond that investigating the actual degree of coherence of 
random lasing emission is gaining weight now HUH]. Besides 
coherence we discuss the scattering mean free path l s of ran- 
dom lasers. In literature that scale has often been estimated to 
be a material characteristics which was independent of the po- 
sition in the sample. We show that this is actually not the case. 
The scattering mean free path itself is intrinsically tied to self- 
consistent nonlinear gain and therefore the spatially dependent 
dissipative character iH, 12] of the random laser, especially at 
the surface. 




FIG. 1: ZnO spherical Mie scatterers at random locations (blue) 
are optically pumped from above (wide yellow beam). The pumping 
yields an inversion of the atomic occupation number within the ZnO 
causing stimulated emission of light (orange light paths). The emit- 
ted intensity multiply scatters and concentrates due to the samples 
density distributions. At the laser threshold the system experiences 
a phase transition and second order coherent intensity escapes the 
system through its surfaces (orange cone). The scatterers radius is 
ro = 600nm, A = 723nm and the samples finite dimension is of 
d = 32fiin. 



MODEL 



The system we consider consists of a randomly scattering 
medium in the form of the simplest slab geometry |7]. This 
slab is finite d sized in the z-dimension and assumed of to be 
of infinite extension in the (x, y) plane. In experimentally rel- 
evant situations this refers to film structures of thickness up 
to 300 prni and an in-plane extension of many wavelengths of 
the used light, i.e. infinitely large. The spherical Mie scat- 
terers J2I 0D are embedded in a homogeneous host material 
which is in presence considered to be passive. Both, scat- 
terer and host medium, are described by means of a com- 
plex dielectric function e s and e?,, respectively. The scatter- 
ers are modeled to be optically active ZnO, with the refrac- 
tive index to be n = 9.1. This sample is optically pumped 
in order to achieve a sufficient electronic population inversion 
within the active medium of the scatterers by means of an inci- 



dent pump laser, perpendicular onto the (x, y)-surface of the 
random laser. The laser feedback is guaranteed by multiple 
scattering. The same mechanism actually supports stimulated 
emission and hence coherent light intensity within the setup. 
The so generated laser intensity then may leave the sample 
through both open surfaces of the sample geometry, the dis- 
sipation channels. The emitted light is eventually observed at 
the surface of the sample in the form of lasing spots which 
comprise to a lasing mode. These lasing modes are of a char- 
acteristic size depending on system parameters such as scat- 
terer size, wavelength of the pump source, filling fraction and 
film thickness etc.. 

This mode size, i.e. the intensity correlation or coherence 
length of the coherently emitted laser intensity, is derived by 
means of field theoretical approach based on the electronic 
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localization theory by Vollhard and Wolfe 



THEORY 

The equation of motion for the electric field of stimulated 
emitted light \& w (r ) within the sample is given by the wave 
equation 



-=-e(r)# u (r) + V 2 ^ 



47T 

= -iu-oj u (r) , (1) 



where we denote c to be the vacuum speed of light and 
jui{r) the external source. The dielectric constant e(r ) = 
€b + Ae V(r), where the dielectric contrast has been defined 
according to Ae = e s — £5. It describes the random arrange- 
ment of scatterers in terms ofV(f) — J2p &r — R ), with 
Sfi (f ) a localized shape function at random locations R. The 
intensity correlation function is then related to the field-field- 
correlation function $, often referred to as the four-point- 
correlation, $ = (*£(r, t)^*(f' ,t' )). Here, the angular 
brackets (...) refer to the disorder average or ensemble aver- 
age of this random system. In order to calculate the field-field- 
correlation <I> the Green's function formalism is best suited. 

The (single-particle) Green's function is related to the 
(scalar) electrical field, Eq. (Q]) by 



dV / M' G{?,?'\t,i! ti). (2) 



In order to study the transport of the above introduced field- 
field-correlation we consider the 4-point correlation function, 
defined now in terms of the non-averaged Green's functions, 
i.e. the retarded G and the advanced Green's function G*, 
where now we find $ ~ (GG*). 

The intensity correlation obeys an equation of motion itself, 
the Bethe-Salpeter equation (BS), given in coordinate space 
given as 



center-of-mass momentum Q|| = (Q x , Q y ) variables. How- 
ever, the finite z-coordinate of the slab is transformed into rel- 
ative and center-of-mass real-space coordinates, i.e. z and Z 
respectively. In this representation only the relative coordinate 
is Fourier transformed. This procedure is justified because the 
relative coordinates of the intensity correlation are related to 
the scale of the oscillating electric field, whereas the center- 
of-mass coordinates are related to the scale of the collective 
behavior of the intensity, which is a significantly larger scale. 
Given that the thickness of the slab is much larger than the 
wavelength of the laser light as discussed above, a Fourier 
transform with respect to this perpendicular relative coordi- 
nate is perfectly acceptable. 

In this representation the BS equation, Eq. (01, may be 
rewritten according to 



AS + 2Reewfi - Aew 2 -2p\\-Q\\+ 2ip z d z 
x^(Z,Z')=AG p (Q ll ;Z,Z , )5(p- P ') + (4) 

E AG ^n) JjSf^M z > z ') 

where we used the abbreviation AG = G — G*. The 
rewritten BS equation, Eq. (0), also known as kinetic equa- 
tion, therefore is seen to be a differential equation in finite 
center-of-mass coordinate Z along the finite thickness of the 
slab. This differential equation is to accompanied by respec- 
tive boundary conditions accounting for the reflectivity of the 
sample surfaces. Eq. (0| is solved in terms of an expansion of 
the correlation $ into its moments, identified as energy den- 
sity and energy current density correlation, respectively. A 
self-consistent expression for the diffusion constant is derived, 
accompanied by a pole structure within the energy density ex- 
pression. 

The lasing behavior in terms of the atomic occupation num- 
bers is incorporated by means of the following laser rate equa- 
tions 



(3) 



§(n,ri ; r 2 ,r 2 )=G R (r 1 ,r , 1 )G A (r 2 ,r 2 ) 

+ G R (r u r 5 )G A (r 2 ,r 6 ) x 

'Xj{r5,r3;r 6 ,r 4 )$(r 3 ,r' 1 ;r4,r 2 ). 



In the BS, we introduced the irreducible vertex function 
l( r 5i r 3i r 6, t a) which represents the scattering interactions 
of the intensity correlation inside the disordered medium of 
finite size Jit]. The irreducible vertex is discussed in the given 
reference in detail. 

To account for the particular form of the system geometry, 
Wigner coordinates are chosen, in which a full Fourier trans- 
form of the spatial coordinates within the infinitely extended 
(x, y)-plane is used, where we use relative gji = (q x ,q y ) and 
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where N = Ni(r,t), i = 0, 1, 2, 3 are the population 
number densities of the corresponding electron level; N to t is 
the total number of electrons participating in the lasing pro- 
cess, jij = 1/tjj are the transition rates from level i to j, 
and 7„ r is the non-radiative decay rate of the laser level 2. 
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7p = 1/rp is the transition rate due to homogeneous, con- 
stant, external pumping. Further n p h = N p h/N to t is the pho- 
ton number density, normalized to N tot . 

We are interested in the stationary limit, i.e. d t Ni = 0, 
hence the above system of equations can be solved for the 
population inversion n 2 = N2/N to t to yield 



n 2 = 



7P + Inr + 721 (n p h + 1) ' 



(10) 



where it was assumed that 732 and 710 are large compared 
to any other decay rate. 

In a last step, the laser rate equations are coupled to the mi- 
croscopic transport theory by identifying the growth term in 
a photon diffusion with corresponding growth term in the de- 
rived equation for the energy density correlation. The equality 
is ensured by finding an approbate imaginary part of the scat- 
terers dielectric function e s for any particular light frequency 
uj and any position Z within the slab geometry. 
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FIG. 2: Calculated coherence length £ of the random lasing modes 
as a function of the calculated scattering mean free path l s . Both 
length scales are given in units of the scatterer radius rn. Differ- 
ent curves correspond to different strengths of the pump intensity P, 
given in units of transition rate 721. The different points along a 
given curve correspond to different filling fractions of the zinc oxide 
scatterers. The points from left to right correspond to filling fractions 
of 60%, 50%, 45%, 40%, and 35%. 



DISCUSSION 

The here developed theory of random lasing includes the 
regular self-consistency of the Vollhard- Wolfe type for the 
diffusion coefficient including the interference effects of emit- 
ted light intensity. The self-consistent results for the imagi- 
nary part of the dielectric coefficient of the laser active scatter- 
ers are derived by coupling the mesoscopic transport to semi- 
classical laser rate equations. The single particle self-energy 
E(w) entering the single particle Green's function, as in Eq. 
(Q3, is approximated as the volume filling fraction of scatter- 
ers in the host medium multipled with the scattering matrix of 




FIG. 3: Calculated scattering mean free path l a across the depth Z of 
the random lasing film for various filling fractions as indicated in the 
legends. The pump rate for all systems is set to be P — O.5721. The 
dependence of l s on the depth originates basically from the depth 
dependent atomic inversion, i.e. on the depth dependent imaginary 
part of the dielectric function e s , which has been self-consistently 
calculated. 



the Mie scatterer. The internal resonances establish threshold 
friendly parameters. 

In Fig. [2] the calculated correlation length £, given in units 
of the scatterer radius r , is presented as a function of the scat- 
tering mean free path l s of the random system. The bullets on 
the curves correspond to several filling fractions. We find, that 
the coherence length increases with the mean free paths and 
again decreases with higher pump intensities. That behavior 
can be interpreted as self-balancing of sample energy at the 
threshold. Also in Fig. [2] Z s is given in units of the scatterer ra- 
dius ro-The different curves correspond to different strengths 
of the pump intensity P, given in units of transition rate 721. 

The displayed points along one graph correspond to numer- 
ical evaluations for different filling fractions of the zinc oxide 
scatterers. In particular, the points from left to right corre- 
spond to filling fractions of 60%, 50%, 45%, 40%, and 35%. 
The correlation or coherence length is found to decrease with 
increasing pump intensity. Further, the dependence on the 
scattering mean free path is decreased for stronger pumping. 
This is in agreement with recent experimental results, as dis- 
cussed in reference Jlj]. 

The discussed coherence length at the surface of the sam- 
ple, Fig. |2] has also been calculated across the sample thick- 
ness. In Fig. |U the coherence length £ for a volume filling 
fraction of the scatterer of 40% and a pump rate of P = O.I71 
is shown in the upper panel. This is accompanied by the cal- 
culated photon number density in the middle panel as well as 
the scattering mean free path l s as a function of the sample 
depth Z. The variation of the mean free path is especially 
noteworthy since it is solely dependent on the different mate- 
rial inversion in center of the film and close to the surfaces. 
Depth dependent measurements of the scattering mean free 
path of pumped laser active media have not been reported so 
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FIG. 4: Calculated coherence length { , photon number density n p h 
and scattering mean free path l s across the slab geometry from sur- 
face to surface. The parameter set is a filling fraction of 40% and a 
pump rate of P — O.I721. 

far. 

CONCLUSION 

In summary, we considered a system of randomly posi- 
tioned laser active ZnO scatterers and we developed a sys- 
tematic theory for random lasing in finite sized disordered 
systems. The self-consistent theory of light localization, in- 
cluding interference effects of the emitted laser intensity, is 
coupled to the laser rate equation in order to obtain the corre- 
sponding nonlinear gain of the lasing system. The calculated 



averaged correlation or coherence length of the occurring las- 
ing spots is systematically studied and found to exhibit a char- 
acteristic dependence on the scattering mean free path of the 
sample in dependency to the in-depth position. These results 
are in agreement with recent experimental findings. 
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